Transverse Phase Locking for Vortex Motion in Square and Triangular Pinning Arrays 
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We analyze transverse phase locking for vortex motion in a superconductor with a longitudinal 
DC drive and a transverse AC drive. For both square and triangular arrays we observe a variety 
of fractional phase locking steps in the velocity versus DC drive which correspond to stable vortex 
orbits. The locking steps are more pronounced for the triangular arrays which is due to the fact 
that the vortex motion has a periodic transverse velocity component even for zero transverse AC 
drive. All the steps increase monotonically in width with AC amplitude. We confirm that the width 
of some fractional steps in the square arrays scales as the square of the AC driving amplitude. In 
addition we demonstrate scaling in the velocity versus applied DC driving curves at depinning and 
on the main step, similar to that seen for phase locking in charge-density wave systems. The phase 
locking steps are most prominent for commensurate vortex fillings where the interstitial vortices 
form symmetrical ground states. For increasing temperature, the fractional steps are washed out 
very quickly, while the main step gains a linear component and disappears at melting. For triangular 
pinning arrays we again observe transverse phase locking, with the main and several of the fractional 
step widths scaling linearly with AC amplitude. 

PACS: 74.60.Ge, 74.60.Jg 



I. INTRODUCTION 

The phenomena of phase locking can occur in systems 
interacting with a periodic potential when an external 
AC drive is superimposed over a DC drive. A resonance 
or frequency matching can occur between the AC fre- 
quency and the internal frequency generated by motion 
over the periodic substrate. The best known example is 
that of Shapiro steps in AC/DC-driven single Joaephson- 
junctionscra and arrays of Josephson-junctionaJ, where 
steps are observed in the current-voltage characteris- 
tics. Shapiro like phase locking has also been, exten- 
sively studied] in charge-density-wave systemsu, spin- 
density wavesB ancL vortices in superconductors with pe- 
riodic substratesU~u. In addition, in elastic media moving 
over random disorder, a washboard signal can be gener- 
ated by the periodicity of the elastic media itselftjO, 
allowing Shapiro like steps to be observed when an AC 
drive is .a]j|p]ied, such as in vortex systems with random 
disordeiOll3. In all these systems the AC drive is su- 
perimposed in the same direction as the DC drive. An- 
other characteristic of Shapiro like phase locking steps 
is that the step width oscillates with AC amplitude and 
frequency between zero and a finite value as an nth order 
Bessel function J„(A), where n = corresponds to the 
depinning or critical current and n > corresponds to 
higher order steps. 

Recently a different kind of phase locking was proposed 
for a 2D system of vortices moving through a periodic 
substrate. Here the AC drive is applied in the perpen- 
dicular disection, transverse to the longitudinal applied 
DC drivetJ. Phase locking steps, distinct from Shapiro 
steps, were observed in the velocity versus applied DC 



drive. The most pronounced step occurs for DC drive 
values that allow a sinusoidal orbit to fit into each square 
pinning plaquette. For larger AC amplitudes a number 
of additional fractional steps were also observed. In addi- 
tion, unlike Shapiro steps, the width of the depinning and 
higher order steps were shown analytically and in simu- 
lation to increase monotonically with AC amplitude as 
{A/uj)'^, where A is the AC amplitude and co is the AC 
frequency. A similar transverse phase locking phenom- 
ena can also occur for an ordered vortex lattice moving 
over random disorder, as was recently demonstratedli3. 
In the random disorder case, the steps differ from those 
observed for square periodic pinning. The step width for 
random disorder is more like a Shapiro step, going as 
oc |Ji(A)|. In addition, for low AC amplitudes, the step 
width grows linearly, rather than quadratically as in the 
case of square pinning. 

There are several open questions for the transverse 
phase locking in the square pinning arrays. It is not 
known how the additional fractional steps scale with AC 
amplitude, what the vortex orbits look like along these 
steps, or what the noise signatures are. It also not known 
if there is any scaling at the transition to the phase locked 
state similar to that found in charge density wave sys- 
tems, where critical scaling is known to occur at depin- 
ning and along-jbha transitions into and out of the phase 
locked regionaljila. The effect of disorder on the trans- 
verse phase locking has also not been investigated. The 
vortex filling fraction or applied field in experiments will 
likely determine the step width or whether the steps oc- 
cur at all. For certain filling fractions the. overall vor- 
tex lattice can be disordered or frustratedo, which may 
destroy the steps. In addition typical transport experi- 
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ments are often done near Tc so the effects of finite tem- 
perature on the fractional and main steps should also 
be investigated to determine the feasibility of observing 
these states experimentally. 

It would also be interesting to examine the effects of 
a transverse AC drive on systems in which there is an 
internal transverse oscillation even in the absence of an 
AC drive. An example of such a system is vortex motion 
through a triangular pinning array. Interstitial vortices 
driven through such an array by a longitudinal DC drive 
oscillate back and forth in the transverse direction as 
they move in order to pass through each minima of the 
interstitial potential. In this case the additional internal 
transverse oscillation may lock with the AC transverse 
force, leading to an enhancement of the transverse phase 
locking similar to that observed for the square case. Al- 
ternatively, the phase locking may be more Shapiro like, 
such as the steps found in the transverse phase locking 
with random pinningEj. 

We focus on the transverse phase locking for vor- 
tex motion in square and triangular pinning arrays in 
thin film superconductors. Vortex pinning and dynam- 
ics in thin samples with periodic pinning have attracted 
considerable attention, since the effects of pinning can 
be systematically controlled. These pitming arrays can 
be Gseated using an array of holescJ^EJ or magnetic 
dota23'c3. Various pipmiig-geometries have also been cce.- 
ated such as-^quareOESo, triangulares, rectangulaitB, 
and KagomeEZI. Pronounced commensurability effects ap- 
pear as peaks in the critical current at integer matching 
fields and fractional fields where the vortices can form a 
symmetrical configuration with the pinning array so that 
the vortex-vortex interaction energy is reduced. These 
matching configurations have been imaged directly in 
experimentE3'LJ and also observed in simulationsEa. For 
large pinning sites above the first matching field, where 
there are more vortices than pinning sites, multiple vor- 
tices can be captured by each pin up to a saturation 
point. Beyond this saturation number, additional vor- 
tices will be located between the pinning sites in the 
interstitial regions. For small pinning sites, the satu- 
ration number is one, so that for all fields above the 
first matching field the vortices are located at the in- 
terstitial regions. Direct imaging and simulations have 
shown that these interstitial. vpr.tices can form ordered 
commensurate arrangementstJcJ. The interstitial vor- 
tices are still pinned even though they do not interact 
directly with the pinning sites, since a periodic poten- 
tial is created by the repulsive interaction of the vortices 
trapped at the pinning sites. The interstitial vortices 
are typically far less strongly pinned than the vortices 
at the pinning sites, so that under an applied drive the 
interstitial vortices move while the vortices at the pin- 
ning sites remain immobile. Evidence for the motion 
of the interstitial yartices has been provided by trans- 
port measurementso as well as direct imaging in square 
pinning arrays, which show that the interstitial vortices 
move in ID paths between the pinning sitescil. In simu- 



lations, ID interstitial vcujt^x flow between pinnina, sites 
has also been observedE3o. Recent experimentsH and 
simulationga have demonstrated the appearance of steps 
in the V-I characteristics for moving interstitial vortices 
when a DC and AC drive are superimposed. In the previ- 
ous theoretical work for an applied transverse AC drive, 
vortex filling fractions between 1.0 and 2.0 were exarg- 
ined, where each pinning site captures only one vortexH. 

In this work we examine the additional fractional steps 
for the square pinning array with a transverse AC drive. 
We find fractional steps both above and below the main 
step. We show that the width of most of the fractional 
steps scales as the square of the AC amplitude, as pre- 
viousljy shown for the depinning current and the main 
stcpO. Along the steps, the vortex orbits are stabilized 
in fixed trajectories, while in the non-step regions, the 
vortex motion is chaotic in appearance, producing a dis- 
tinct velocity noise spectra. We also find scaling in the 
velocity versus applied DC drive at depinning and near 
the main step, with (3 — 1/2. The steps are most promi- 
nent for commensurate vortex filling fractions where the 
interstitial vortices form a symmetrical ground state. For 
filling fractions where the vortices are disordered, the 
phase locking phenomena is absent, and for filling frac- 
tions just off of a commensurate filling, the steps have 
a linear increase with drive due to the fact that certain 
portions of the sample are not undergoing phase locking. 
Wc find that even for small finite temperatures, the frac- 
tional steps become completely washed out. The main 
step is visible up to the melting point of the interstitial 
vortex lattice. With triangular pinning we find a much 
more pronounced main step and fractional steps, with 
the depinning threshold nearly the same as that of the 
square case. The widths of the main and most of the 
fractional steps scale linearly with the AC amplitude for 
triangular pinning. 

II. SIMULATION 

We numerically integrate the overdamped equations of 
motion, which for a single vortex i is 

f, = f™ + + foe + fAC = V, (1) 

The total force acting on vortex i is fi, f^"" is the total 
force from all the other vortices, is the force from the 
substrate, ioc is the force from the DC drive and iAc is 
the force from the AC drive. The simulation is in 2D with 
periodic boundary conditions in the x and y directions. 
For the square pinning case we place an array of Np = 
N X N pinning sites in the sample, where each pinning 
site is an attractive parabolic potential of radius rp and a 
maximum force fp. To create a triangular pinning array, 
we displace every other row of pinning sites by half a 
pinning lattice constant. The force from the pinning sites 
is /r = E^(/pAp)|r. - ri^^l e(r, - |r, - ri^)|)f(^\ 
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DC Drive 

FIG. 1. Velocity V vs applied DC drive fnc for AC am- 
plitudes A = 0.25 (upper curve) and 0.35 (lower curve) for 
a square pinning array. The upper curve has been offset up- 
wards by 0.02 for clarity. 

(v) 

where Q is the step function, r^T is the location of pin- 
ning site fc, and f|^'' = (r — r^'^')/|ri — r^^'']. In our system 
each pin captures only one vortex. The total force from 
the vortex- vortex interaction is f™ = — X^j^i Vi[/u(r.y ) 
where we use U{r) — ~ln{r), the Pearl vortex- vortex 
interaction potential, appropriate for thin film supercon- 
ductors. The experiments with periodic pinning arrays 
are in this 2D limit. FoMomputational efficiency we use 
a summation techniqua23 to evaluate the vortex-vortex 
interaction. Temperature can also be applied to the sys- 
tem by adding a noise term ff to the equation of mo- 
tion. This noise term has the properties < f^{t) 0.0 
and < ff{t)fj{t) >= 2rikBT5^j5{t - t). Here we set 
= fcs = 1. 

We have considered pinning arrays from 4 x 4 to 10 x 
10 pins and observe the same results in each case. The 
initial vortex configurations are obtained by annealing 
from a high temperature where the vortices are in a liquid 
state and cooling to T = 0. The DC driving is in the 
x-direction, ioc = Jdc^, and the AC force is in the y- 
direction, f^c = Asm{ujt)y. The DC force is increased 
from zero in increments of 0.0009. We average the vortex 
velocities = ^ • x at each increment for 175000 MD 
steps, and the resulting DC force versus velocity curve is 
proportional to the DC voltage-current curve. 



III. FRACTIONAL STEPS FOR SQUARE 
PINNING ARRAYS 

In Fig. 1 we show the vortex velocity V versus applied 
DC drive foe for two different AC amplitudes, A = 0.25 
(upper curve) and A = 0.35 (lower curve), for a fixed fre- 
quency oi lli — 0.002. Here there is a large step, centered 
near foe = 0.32 for A = 0.25 and foe = 0.35 




0.222 0.232 0.4 0.42 0.44 0.46 



'do DC 

FIG. 2. Magnified regions of vortex velocity V versus ap- 
plied DC drive foe for A = 0.35. (a) 0.222 < foe < 0.241, 
showing phase locking steps below the main step. (b) 
0.4 < foe < 0.45, showing phase locking steps above the 
main step. 

for A — 0.35, which corresponds to a phase where the 
vortex orbit is commensurate with the pinning, forming 
a sinusoidal like orbit as studied previously. We label 
this step the "main step" or 1:1 step. In addition to this 
step, a number of smaller or fractional steps are visible 
both above and below the main step. For example, for 
A = 0.35 a prominent step appears near fnc — 0.27, and 
a small step is visible near /^c = 0.46. The widths of all 
the steps, including the depinning threshold, are smaller 
for the lower AC amplitude. The depinning threshold 
drops from foe = 0.196 for A = 0.35 to foe = 0.190 
for A = 0.25, while the width A of the main step drops 
from A = 0.09 for A = 0.35 to A = 0.04 for A = 0.25. In 
Fig. 2(a) we show a magnified region of 0.222 < foe < 
0.241, near the depinning threshold, for the A = 0.35 
curve in Fig. 1. Here the large number of fractional steps 
present are clearly visible, with the step widths getting 
progressively smaller for lower foe- In Fig- 2(b) we show 
a magnified region of the A = 0.35 curve above the main 
step, from 0.4 < fuc < 0.45, where a series of fractional 
steps appear. Again the step widths decrease as fuc is 
lowered toward the main step. 



A. Commensurate Vortex Orbits and Noise 

On steps at small foe values, a vortex spends more 
than one AC cycle in a single plaquette, while on steps 
at larger foCi the vortex moves through more than one 
plaquette in a single AC cycle. In Fig. 3 we illustrate the 
vortex trajectories (lines) and the vortex positions (dots) 
for fixed fuc drives on and near the steps in the V versus 
foe curve of Fig. 1 for A = 0.35. Fig. 3(a) shows the 
vortex trajectory for the step near fjjc = 0.235, which 
is most clearly visible in Fig. 2(a). Here, a stable vortex 
orbit occurs in which the vortex spends two AC cycles in 
one plaquette before moving over to the next plaquette to 
repeat the sequence. We observe similar behavior at the 
smaller steps for fuc < 0.235, where in the stable orbits, 
the vortex spends n or n -I- 1/2 complete AC cycles in one 
plaquette. In Fig. 3(b), in the stable orbit for the step at 
foe — 0.243 in Fig. 1, the vortex spends 3/2 of an AC 
cycle in each plaquette. In Fig. 3(c) we show the vortex 
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(e) (f) 
FIG. 3. Vortex trajectories (lines) and vortex positions 

(black dots) over equal time intervals for constant foe values 

chosen from the lower curve in Fig. 1 with A — 0.35. The 

vortices in the square lattice are trapped in pinning sites and 

remain immobile; only the interstitial vortices are moving. 

The trajectory of only a single mobile vortex is presented for 

clarity, foe = (a) 0.235, (b) 0.243, (c) 0.28, (d) 0.305, (e) 

0.36, and (f) 0.429. 

trajectory for fjjc — 0.28, which corresponds to the sec- 
ond largest step in Fig. 1. Here the vortex goes through 
one AC cycle per plaquette. In the non-step region of 
Fig. 3(d), at foe = 0.305, there is no single stable vor- 
tex orbit. Instead the vortex trajectory wanders over 
time. There is an area close to the occupied pinning sites 
where the vortex never flows due to the vortex-vortex 
repulsion. We find similar orbits for the other non-step 
regions. On the main step, illustrated in Fig. 3(e) at 
foe — 0.36, the vortex moves through two plaquettes in 
a single AC cycle. In Fig. 3(f) the vortex orbit for a step 
above the main step, foe = 0.429, has a similar pattern 
to that seen on the main step, Fig. 3(e), except that the 
vortex moves through three plaquettes in each complete 
AC cycle. 

The presence of stable orbits can also be inferred 
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FIG. 4. Velocity noise S{f) versus frequency / for 
A — 0.35. Lower curve: Noise signature along the main step 
corresponding to the vortex orbit at foe = 0.36 shown in 
Fig. 3(e). Upper curve: Noise signature in the non-step re- 
gion for the vortex orbit at foe ~ 0.305 shown in Fig. 3(d). 




FIG. 5. (a) Log-log plot just above depinning of 

V = {f — fc) for A = 0.35, where fc = 0.196 is the depinning 
threshold and / is the DC driving force. The solid line shows 
a slope of /9 = 1/2. (b) Log-log plot above the main step of 

V — Vq — {f — fq) where Vq = 0.24 is the velocity along the 
main step and fq = 0.395 is the DC driving force where the 
step ends. The solid line shows a slope of (3 — 1/2. 

through measurement of the velocity noise, which corre- 
sponds experimentally to the voltage noise. We acquire 
a time series of the velocity in the a;-direction for a fixed 
foe value. The corresponding noise spectrum is given 
by 5'(/) = \ JV{t)e-''^''f*dt\'^. In Fig. 4, we examine the 
noise spectra for A = 0.35 for a step region, foe = 0.36, 
and a non-step region, foe — 0.305, corresponding to 
the orbits in Fig. 3(e) and Fig. 3(d). In the step region 
(lower curve), a strong periodic signal is present at the 
frequency induced by the AC driving, as indicated by 
the narrow large spike and the higher harmonics. For 
the non-step region (upper curve) there is is a broadened 
peak near the AC driving frequency. The noise power 
for the lower frequencies is over 4 orders of magnitude 
higher than for the step region. This lower frequency 
noise power comes from the vortex orbits drifting over a 
time scale that is longer than the AC driving time scale. 
The noise spectra does not have a broad band or 1//" 
characteristic. We find a similar noise spectra for the 
other non step regions. 

B. Scaling Velocity vs Drive and Fractional Steps 

The velocity curves in Fig. 1 show considerable cur- 
vature near the depinning threshold and just above each 
step. Fisher has suggested that depimung can be consid- 
ered a dynamical critical phenomenally where 

V={f- fcf, (2) 

where / is the DC driving force, fc is the driving force 
at which depinning occurs, and V is the velocity of the 
driven media. In addition, in numerical work on CDW 
systejn.s a similar scaling is found near the mode locked 
stepstZEa in the form 

\V-V,\^\f-f,f, (3) 

where Vq is the velocity along the qth step and fq is the 
driving force at which the step ends or begins. For a sin- 
gle particle moving in a ID periodic substrate, P = 1/2. 
In charge density wave systems with random pinning, 
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FIG. 6. The dV/dfoc curve for A = 0.35 starting at the 
DC depinning threshold of foe = 0.196. 

simulations of the depinning and step transitions give (3 = 
0.67 in 2D and 0.83 for 3D, while experimentalxesults are 
consistent with the results for the 3D regimeli3. In our 
system, one might expect to find scaling with (3 = 1/2 as 
in the ID periodic substrate picture since the vortices are 
commensurate with the periodic pinning. On the other 
hand, the orbits in Fig. 3 indicate that the vortices have 
a 2D velocity component and do not flow in strictly ID 
channels, so it is not clear if or how the velocities will 
scale. 

In Fig. 5(a) we plot on a log scale V — f — fc just above 
depinning for A — 0.25, with the depinning threshold of 
/e = 0.196, and in Fig. 5(b) we plot \V - Vg\ = f - fg 
above the main step for A = 0.25, with the velocity Vq = 
0.24 along the main step, and where the end of the step 
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FIG. 7. The scaling of the step width A as a function of 
AC amplitude A for the fractional steps corresponding to; 
(diamonds) foe = 0.235 [Fig. 3(a)], (triangles) foe = 0.243 
[Fig. 3(b)], and (squares) fnc = 0.28 [Fig. 3(cj]. The scal- 
ing of the main step at foe ~ 0.36 [Fig. 3(e)] is also shown 
(circles). The solid line is a quadratic fit. 
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FIG. 8. The V vs fnc curves for different vortex filling 
fractions. The curves have been shifted up from V = 0.0 for 
presentation. From lowest to highest: B/B^ = 17/16, 5/4, 
0.34, 3/2, 0.67, 1.97, 2.0, and 2.04. Here A = 0.35. 

is at fq — 0.395. Here in both cases for low / the curves 
fit well to /3 = 1/2, with a deviation to linear behavior for 
higher /. Our velocity resolution was too low to obtain 
adequate scaling for steps other than the main step. We 
find similar scaling for the depinning and the main step 
for all filling fractions at which the interstitial vortices 
form a symmetrical configuration. For non-symmetrical 
vortex configurations, the scaling breaks down due to the 
fact that the depinning process for the non-symmetrical 
configurations is plastic, rather than elastic as in the case 
of the symmetrical configurations. Scaling may still occur 
for the incommensurate systems as found in simulations 
with random pinning where plastic flow occurs j-where 
scaling in the velocity versus drive was observedEll with 
/3 > 1.5. Very large systems beyond the scope of this 
work would be needed in order to determine if such scal- 
ing occurs for the incommensurate cases. In addition, 
some of the scaling breaks down for high AC amplitudes, 
as seen for A — 0.35 in Fig. 1 where there are sharp jumps 
into and out of certain steps. For example, in Fig. 6 we 
plot the dV/dfuc curve for A = 0.35, showing peaks go- 
ing into and out of the phase locked regions. For /3 = 1/2 
scaling, the overall shape of dV/dfoc should scale with a 
power —1/2, which is consistent with the curve in Fig. 6. 

We next show how the widths of different steps scale 
with AC amplitude. In earlier work it was shown that 
the increase in the depinning force and the width of the 
main step increase as the square of AC amplitude. In 
Fig. 7 we plot the step widths A vs AC amplitude A 
corresponding to the steps at Jdc — 0.235 [Fig. 3(a)], 
foe = 0.243 [Fig. 3(b)], and foe = 0.28 [Fig. 3(c)]. 
We also show the scaling for the main step, foe = 0.36 
[Fig. 3(e)], for a comparison. The widths all fit well to 
A^ with some deviations at higher AC amplitudes, where 
the widths increase faster than A^ . 



5 




FIG. 9. Vortex trajectories (lines) and positions (black 
dots) along the main step for fnc = 0.36 and A — 0.35 
for (a) B/B^ = 2.0, showing complete phase locking; (b) 
B/B^ = 2.04, showing that only a certain portion of the 
moving vortices are undergoing phase locking. 



ground state for B / = 17/16, 5/4, 3/2, and 2.0. At 
these fillings the phase locking steps are the most pro- 
nounced, and the velocity curves are almost identical for 
each one. These are also the same filling fractions where 
peaks in the critical current are observed in experiments 
and simulations. Since the interstitial vortices form a 
symmetrical ground state, the effective interstitial vortex 
interactions cancel, and the system can be thought of as 
a single particle moving through a periodic substrate. In 
the absence of an AC drive, the vortices collectively move 
in ID paths under a DC drive. For the filling fractions 
B/B^ — 0.34 and 0.67, the ground state is disordered 
and the phase locking is absent. In the absence of an 
AC drive the vortices move plastically in winding paths. 
Under an applied AC drive, this winding motion destroys 
any phase locking. 

For filling fractions B /B^ = 1.97 and 2.04, just below 
and just above the commensurate filling at B /B^p = 2.0, 
the phase locking still occurs; however, the steps have 
an additional linear velocity increase superimposed on 
them. In Fig. 9(a) we illustrate vortex trajectories along 
the main step for B/B^ = 2.0 and in Fig. 9(b) for 
B/B^ = 2.04. For B/B^ = 2.0 the vortices move 
in well defined stable sinusoidal trajectories, while for 
B/Bfj, = 2.04 only certain rows of vortices have this or- 
dered motion while the vortices in the other rows have 
unstable ordering. Only the rows with stable orbits are 
phase locked and remain at the same velocity along the 
step for increasing foe- The vortices moving through 
the unstable channels are not phase locked and there- 
fore their velocity increases with foCi accounting for the 
linear velocity increase in the V - fnc curves along the 
step. The stable moving vortex rows have a commensu- 
rate number of vortices. In the system presented here, for 
an 8 X 8 pinning lattice there are eight interstitial vortices 
in each commensurate row. There are more than eight 
vortices in the rows with unstable vortex trajectories, so 
the rows are incommensurate. For certain filling frac- 
tions such as B/B^ — 0.67, all the rows are incommen- 
surate and there is no phase locking. For B/B^ = 1.97, 
where the unstable or incommensurate rows have less 
than eight vortices, we observe vortex trajectories simi- 
lar to those in Fig. 9(b). We note that similar behavior 
occurs just above and below the filling fractions of 5/4 
and 3/2. These results are consistent with experimental 
results in which Shapiro steps were most clearly defined 
at B/B^ = 2.d3. 

D. Effect of Temperature 



C. Effect of Filling Fraction 

In Fig. 8 we show a series of velocity vs fnc curves for 
increasing filling fractions from (bottom to top) B/B^ = 
17/16, 5/4, 0.34, 3/2, 0.67, 1.97, 2.0, and 2.04. The 
interstitial vortex lattice forms an ordered symmetrical 



We consider the effects of temperature by adding a 
noise term to the equation of motion for the vortices. 
We normalize our temperature by the melting tempera- 
ture Trm where is the temperature at which the vor- 
tices began to diffuse in the absence of an external drive. 
In our system this onset is sharp and well defined, due 
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FIG. 10. The V vs foe curves for different tempera- 
tures for B/B^ = 2.0 and A — 0.35. The curves have been 
shifted up from zero for presentation. From lowest to highest, 
T/Tm = 0.0, 0.07, 0.28, 0.64, 0.75, 0.87, and 0.96, where T™ 
is the melting temperature. 

to the fact that the vortices are sitting in a periodic sub- 
strate. In the absence of a substrate the onset of diffusion 
is more gradual and a well defined melting temperature 
does not appear. Experimentally, the meltiiig|-©f the in- 
terstitial vortices occurs at or just below TcEili^J so that 
Tm ~ Tc. In Fig. 10 we show a series of V vs foe curves 
for A ^ 0.35 for T/T™ = 0.0, 0.07, 0.28, 0.64, 0.75, 0.87, 
and 0.96, from bottom to top. The fractional steps wash 
out very quickly with temperature at about T/Tm = 0.1. 
The main step is visible all the way up to r,„ but is more 
cusp like for T/r,„ > 0.64. This washing out of the step 
with temperature is also consistent with the Shapiro step 
experiments near Tc which found only a cusp featured. 
Experimentally it would still be possible to measure the 
increase of the step width with AC amplitude by taking 
the derivative of the I ~ V curves, which shows a dip at 
each side of the step. Since most transport experiments 
are performed close to Tc, it will be very difficult to ob- 
serve the fractional steps. The widths of the steps can be 
increased by increasing A, allowing fractional steps to be 
visible for T > QATm- In practice, however, there will be 
a limit to how large A can be made before the vortices 
at the pinning sites begin to depin. 

IV. TRANSVERSE PHASE LOCKING FOR 
TRIANGULAR PINNING 

Next we study the transverse phase locking for the 
case where the pinning is triangular rather than square. 
We use the same system size and pinning density as for 
the square pinning array, but every other pinning row is 
shifted in the a:-direction by half a lattice constant. Here 
the motion of an interstitial vortex for zero AC amplitude 



0.8 



0.6 



>0.4 



0.2 " 








0.2 0.4 

^DC 

FIG. 11. The velocity V vs DC driving force foe curves 
for a triangular pinning array at different AC amplitudes A. 
The system size and pinning density are the same as for the 
system in Fig. 1. The curves have been shifted up for pre- 
sentation. From lowest to highest, A = 0.35, 0.25, 0.15, 0.05, 
and 0.0. 

has a periodic transverse velocity component, unlike the 
square case where for zero AC amplitude the motion is 
strictly ID. In Fig. 11 we show V vs foe for different 
values of A. We find phase locking similar to that seen 
in the square pinning case. The step widths all increase 
monotonically with AC drive, rather than oscillating as 
expected for Shapiro steps. There are some noticeable 
differences from the square case, in that the depinning 
force is only weakly affected by increases in ^. In addi- 
tion the main step is wider than that observed for square 
pinning at the same AC amplitudes. We also find that the 
step widths increase linearly with AC amplitude, rather 
the quadratically as in the square pinning array. 

In Fig. 12(a) we show the vortex motion above depin- 
ning for A = 0.0, indicating that the vortex moves with 
a periodic transverse component even in the absence of 
an AC drive. In Fig. 12(b-f) we show the vortex or- 
bits along various steps for A — 0.25. In Fig. 12(b), 
at foe = 0.07, a stable periodic orbit forms where the 
vortex goes through an interesting loop that alternates 
from up to down. This loop feature appears for stable 
orbits where foe < 0.1. For the step near fjjc = 0.11 
[Fig. 12(c)], corresponding to the second largest step, a 
stable orbit forms similar to that in Fig. 12(b) but the 
loop feature is absent. On the small step at foe = 0.1375 
[Fig. 12(d)], a stable orbit forms which is not symmetri- 
cal in the y-direction. A vortex trajectory for a non-step 
region at foe = 0.2 [Fig. 12(e)] has similar characteris- 
tics as the trajectories for the non-step regions for square 
pinning. In Fig. 12(f) we show the vortex orbit along the 
main step, fuc = 0.3, where a well defined sinusoidal or- 
bit can be seen as the transverse frequency generated by 
the DC motion matches one-to-one with the frequency of 
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(a) (b) 



(c) (d) 



(e) (f) 
FIG. 12. Vortex trajectories for triangular pinning on vari- 
ous steps for (a) A = 0.0 and (b-f) A = 0.25. (a) For A = 0.0 
above depinning, foe = 0.3, there is a periodic transverse 
motion of the vortex even for zero AC drive. Remaining pan- 
els: A = 0.25 and (b) foe = 0.07, (c) foe = 0.11, (d) 
foe = 0.1375, (e) foe = 0.2, and (f) foe = 0.3. 

the transverse AC drive. For higher DC drives, there 
are additional stable orbits along the steps as in the case 
for the square pinning. Similar noise signals for step and 
non-step regions as seen in Fig. 4 are also observed for 
the triangular pinning case. 

In Fig. 13 we show the scaling of the width A of the 
main step for triangular pinning with AC drive A, along 
with the width A of the main step for square pinning for 
comparison. Here the magnitude of A for the triangular 
pinning is greater than that for the square pinning for all 
AC values considered here. The width of the triangular 
steps also scales linearly with AC amplitude rather than 
quadratically as in the square pinning, so that for smaller 
AC amplitudes the difference in the magnitudes of A 
becomes greater. This result suggests that the detection 
of transverse phase locking would be easiest for a system 
with triangular rather than square pinning arrays. 

10°, . . ^^^^ 
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FIG. 13. Scaling of the step width A versus AC amplitude 
A of the main step for triangular pinning (squares), along with 
the scaling for the main step for square pinning (circles) for 
comparison. The solid line is a linear fit. 
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FIG. 14. Scaling of the step width A vs y4 for triangular 
pinning for (circles) foe ~ 0.11 [Fig. 12(c)], along with (tri- 
angles) the depinning force. We also show the scaling of the 
main step (squares) for comparison. 

In Fig. 14 we show the scaling of the step width A with 
AC amplitude A for the fractional step at foe = 0.11 
(the orbit shown in Fig. 12(c)), along with the width of 
the main step and the change in value of the depinning 
force. Here the fractional step scales linearly with AC 
amplitude, just as the main step does. The difference 
in the depinning force at different values of A, which 
is much smaller than the widths of the main step and 
the fractional step, does not scale linearly with A but 
increases quadratically or even faster. 

We can compare our results for transverse phase lock- 
ing in the triangular pinning— to the results for random 
pinning. For random pinningEj the vortices in the highly 
driven phase regain partial triangular ordering and a 
washboard signal appears. Unlike the square pinning 
case, in random pinning in the absence of a transverse 
AC drive the vortices do not move in strictly ID chan- 
nels but have a transverse velocity component centered 
around zero which can also have a washboard signal. For 
the random pinning it was found that the width of the 
transverse phase locking step oscillates as a function of 
AC amplitude, which is different from the results for the 
triangular or square pinning. For low AC amplitudes, 
however, the increase in the step width is linear with AC 
amplitude in the random pinning case. Using a simple 
model it was shown that any form of disorder that in- 
duces a transverse temporal order in the absence of an 
AC drive gives rise to a linear-dependence of the step 
width for small AC amplitudesCj. This is similar to the 
case of the triangular pinning, where there is transverse 
temporal order present in the zero AC drive limit. In 
addition, for the case of random pinning, the vortex tra- 
jectories along the step show stable sinusoidal orbits in a 
similar manner to the stable orbits found in the periodic 
pinning cases. 
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V. CONCLUSION 

We have investigated the transverse phase locking in 
square and triangular pinning arrays. For square arrays 
we find a series of fractional steps which correspond to 
stable vortex orbits. Along the non-step regions there 
are no stable periodic vortex orbits. All the fractional 
steps increase monotonically in width with AC ampli- 
tude. Along the steps a narrow band noise signal is 
present, while in the non-step regions the narrow band 
signal is broadened and large low frequency noise power 
appears. We have shown explicitly that the width of sev- 
eral of the fractional steps increases as the the square of 
the AC amplitude. The velocity vs drive near the depin- 
ning threshold and the main step scale with (3 = 1/2. 

We find that the phase locking is most pronounced for 
interstitial vortex filling fractions at which the intersti- 
tial vortices form a symmetrical ground state, such as 
at B/Bff, = 1.25,1.5 and 2.0. For filling fractions near 
these commensurate fillings, a partial phase locking oc- 
curs where certain regions of the sample have stable phase 
locked orbits while other regions are unstable and the 
vortex velocity in these regions increases linearly. For 
filling fractions where the ground states are disordered 
the phase locking is absent. With a finite temperature, 
the fractional steps appear only at low temperatures and 
are washed out at higher temperatures, while the pres- 
ence of the main step can be detected up to the melting 
temperature Tm- 

For triangular pinning arrays, where moving intersti- 
tial vortices have a periodic transverse motion even in 
the absence of a transverse AC drive, we again observe 
a transverse phase locking with all the step widths in- 
creasing monotonically with AC amplitude rather than 
oscillating as in the case of Shapiro steps. The depinning 
force is only weakly increased by increasing the AC am- 
plitude. The main step for the triangular case is larger 
than that observed in the square pinning case, and its 
width increases linearly with AC amplitude, unlike the 
quadratic increase observed for the square case. The lin- 
ear increase in the step widths with AC amplitude is 
similar to the results for transverse phase locking with 
random disorder with small AC drives, where there is 
also an ordered transverse motion in the absence of an 
AC drive. We show that the width of the fractional steps 
also increase linearly with AC amplitude, but the differ- 
ence in the depinning force increases quadratically. 

Our predictions should be testable for superconductors 
with periodic pinning arrays where only one flux line is 
captured per pinning site. The steps can best be ob- 
served in samples where Shapiro steps have already been 
observed with longitudinal DC and AC drives at filling 
fractions where a symmetrical vortex configuration oc- 
curs, such as B/B^p = 17/16, 5/4, 3/2, or 2. For experi- 
ments performed near Tc, it is unlikely that the fractional 
steps can be observed; however, the main steps should be 
visible. A particular advantage of transverse phase lock- 



ing steps over longitudinal (or Shapiro) steps is that the 
step width can be made arbitrarily large by increasing the 
AC amplitude or lowering the AC frequency. Our results 
should also be relevant for vortex motion in Joscphson- 
junction arrays at commensurate fillings. 
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